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ABSTRACT 


PAFT  177,  a  program  for  the  analysis  of  fault  trees  coded  In  Fortran  77,-4&  - 
presentedr  Given  the  structure  of  a  fault  tree  and  the  probability  or 
failure  rates  of  its  basic  events,  PAFT  F77  calculates  the  probabilities 
of  the  top  and  all  intermediate  events,  as  well  as  the  marginal  Importances 
of  all  basic  events.  When  the  input  is  in  terms  of  failure  rates  Instead 
of  probabilities,  the  marginal  occurrence  rates  of  all  basic  events  and  the 
occurrence  rate  of  the  top  event  are  also  calculated.  This  additional  in¬ 
formation  though,  is  invalid  for  fault  trees  comprising  NOT  gates  or  those 
gate  types  that  must  be  modeled  with  NOT  gates  such  as  exclusive  OR  gates. 
The  program  is  designed  to  handle  simple  fault  trees  and  those  with  repli¬ 
cated  basic  events. 

The  introduction  to  this  text  describes  the  program's  ability  and  limita¬ 
tions  together  with  a  brief  theoretical  review  of  fault  tree  analysis.  The 
second  section  is  a  users  guide  to  PAFT  F77  on  the  UNIX  operating  system  at 
the  University  of  California,  Berkeley.  The  final  section  gives  an  in  depth 
description  of  the  program’s  Fortran  77  code. 
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INTRODUCTION  TO  PAFT  F77 


1.1  What  is  a  Fault  Tree? 

The  safety  of  nuclear  reactors  Is  one  of  the  many  engineering  risk 
considerations  In  the  forefront  of  today's  political,  environmental  and 
social  dialogue.  Answering  these  questions  of  safety  from  an  emotional  or 
political  point  of  view  falls  short  of  the  scientific  objectivity  which  our 
society  expects  and  demands.  Fault  tree  analysis  Is  a  method  available  to 
the  engineer  for  determining  reliability  of  complex  systems.  A  fault  tree 
Is  a  directed  acyclic  graph  representing  the  relationship  between  a  system 
event  and  subsystem  events.  In  turn,  each  subsystem  event  can  be  represented 
by  more  detailed  subsystem  events  until  fundamental  or  basic  events  are 
reached.  The  state  of  each  fundamental  component  or  basic  event  will  de¬ 
termine  the  success  or  failure  of  the  entire  system,  through  the  logic 
structure  of  the  fault  tree. 

Constructing  a  fault  tree  for  a  complex  system  can  be  a  formidable 
task.  The  millions  of  components  of  the  system  must  be  identified  and 
mathematically  modeled.  Then  this  data  must  be  used  to  construct  a  fault 
tree  which  correctly  represents  the  logical  dependence  between  the  system's 
component  events.  For  the  fault  tree  of  a  nuclear  reactor,  the  top  event 
could  be  the  catastrophic  occurrence  of  core  meltdown,  while  a  basic  event 
could  be  as  Insignificant  as  the  outside  humidity.  Fault  tree  analysis  Is 
not  limited  to  nuclear  reactor  reliability  but  Is  a  commonly  used  statis¬ 
tical  and  engineering  tool.  The  following  example  should  give  some  idea 
of  the  variety  of  applications  of  this  method. 

You  are  offered  a  game  of  chance  involving  a  fair  coin.  If  either  of 
the  first  two  flips  Is  a  head  and  either  of  the  second  two  flips  Is  a  head 
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you  lose  $1,  otherwise  you  win  $1.  In  this  example,  the  top  event  is  not  a 
catastrophe,  but  still  the  basic  principles  of  probability  behind  this  game 
of  chance  can  be  applied  to  more  serious  and  consequential  cases ,  The  fault 
tree  (in  this  case  It  would  more  likely  be  called  a  game  tree)  Is: 


Lose  $1 
AND  gate 


Heads  Heads  Heads  Heads 

Basic  event  //I  Basic  event  #2  Basic  event  #3  Basic  event  //4 

The  probability  of  each  occurrence  of  a  head  is  The  probability  of 

3  9 

each  OR  gate  is  and  the  probability  of  the  AND  gate  or  top  event  is 

Therefore,  with  the  opportunity  quoted,  one  would  gracefully  decline  this 

game  since  your  expected  winnings  are  --g  of  a  dollar. 

1.2  The  Program's  Ability 

The  input  to  PAFT  F77  is  the  structure  of  a  fault  tree  together  with 
the  probabilities  or  failure  rates  of  its  basic  events.  The  program  returns 
the  probabilities  of  all  gate  events  including  the  top  event  and  the  marginal 
liiq>ortances  of  all  basic  events.  The  marginal  importance  of  a  basic  event 
is  the  partial  derivative  of  the  probability  of  the  top  event  with  respect 
to  the  probability  of  the  basic  event.  The  marginal  importance  can  also  be 
interpreted  as  the  probability  of  the  top  event  given  that  the  basic  event 
has  occurred,  minus  the  probability  of  the  top  event  given  the  basic  event 
has  not  occurred.  When  analysing  a  system,  it  is  useful  to  know  which  basic 
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events  most  Influence  the  occurrence  of  the  top  event.  The  marginal  impor¬ 
tance  is  one  measure  of  such  Influence. 

When  the  input  is  in  terms  of  failure  rates  Instead  of  probabilities, 
the  marginal  occurrence  rates  of  all  basic  events  and  the  occurrence  rate 
of  the  top  event  are  also  calculated.  The  marginal  occurrence  rate  of  a 
basic  event,  like  its  marginal  Importance,  is  a  measure  of  the  Influence 
a  basic  event  has  on  the  top  event.  The  occurrence  rate  of  the  top  event 
IS  the  sum  of  all  of  the  marginal  occurrence  rates  of  the  basic  events. 

When  a  fault  tree  contains  NOT  gates  or  those  gate  types  that  must  be 
modeled  with  NOT  gates  such  as  exclusive  OR  gates,  the  program,  which  as¬ 
sumes  monotonlclty  of  the  top  event  with  respect  to  basic  events,  calculates 
the  marginal  occurrence  rates  of  the  basic  events  under  these  gates  in¬ 
correctly.  Furthermore,  the  stim  of  these  quantities  produces  an  incorrect 
occurrence  rate  for  the  top  event.  Hence,  marginal  occurrence  rates  may  be 
of  little  value  for  fault  trees  with  NOT  gates  and  exclusive  OR  gates. 

1.3  Theoretical  Background 

The  following  functions  are  used  to  calculate  the  probabilities  of  AND, 
OR,  EXOR,  and  NOT  gates. 

^’and  “  ”  Pi 

PoR  -i-n(i-p^) 

^EXOR  “  Pi  Pz  ■  ^PlPz 

^’not  -  1  -  p  • 


The  remaining  gate  types  (such  as  two  out  of  three  (2/3)  or  exactly  5  out  of 
7  (x5/7))  use  the  following  generating  function. 
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n 

n  (q.  +  p.z)  . 

i-1 

2  3 

For  the  two  out  of  three  case,  the  sum  of  the  coefficients  of  z  and  z 
is  the  correct  probability  with  n  *  3  ;  for  the  event  exactly  5  out  of  7, 
the  coefficient  of  z^  is  the  correct  probability  with  n  ■  7  . 

The  marginal  importance  of  a  basic  event  is 

P(top  I  basic  event)  -  P(top  |  not  basic  event)  . 

The  marginal  occurrence  rate  is  (at  a  specified  time  t) 

(failure  rate  of  basic  event) (marginal  Importance  of  basic  event) 

(1  -  P(baslc  event)) /(I  -  P(top))  . 

If  r^(t)  is  the  occurrence  rate  of  basic  event  i  at  time  t  , 
F^(t)  is  the  probability  event  i  occurs  (for  the  first  time)  by  time 
t  ,  Ij^(t)  is  the  marginal  Importance  of  basic  event  i  at  time  t  and 
P(T  >  t)  is  the  probability  the  top  event  has  not  occurred  by  time  t  , 
then  the  basic  event  occurence  rate  at  time  t  is 

rj^(t)[l  -  F^(t)3l^(t)/P(T  >  t)  . 

If  input  failure  rates  are  constant  in  time,  then  the  distribution  of  time 
to  first  occurrence  of  a  basic  event  is  exponential. 

The  occurrence  rate  of  the  top  event  is 

all 

basic  events 

^  marginal  occurrence  rate  . 

Again  the  formula  for  marginal  occurrence  rates  holds  only  for  trees  not 
comprising  NOT,  EXOR,  or  exclusive  x  out  of  y  gates. 
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1.4  The  Program's  Limitations 

PAFT  F77  Is  an  exponential  time  process  In  terms  of  the  number  of 
replicated  basic  events  contained  In  the  Inputed  fault  tree.  Thus,  the 
program  Is  limited  to  trees  with  a  few  replicated  basic  events,  approxi¬ 
mately  ten  or  less.  The  maximum  dimension  of  an  inputed  tree  is  ninety-nine 
gates,  including  all  basic  events.  Of  these  ninety-nine,  a  maximum  of  20 
replicated  basic  events  can  be  specified,  though  a  tree  with  this  number 
will  take  days  to  compute.  Replicated  gates  are  not  allowed. 

A  second  limitation  to  the  program  Is  Its  Inability  to  handle  NO 
EXOR,  or  exclusive  x  out  of  y  hates  properly  when  calculating  mar{  *'1 
occurrence  rates.  An  enhancement  to  solve  this  problem  can  be  readll- 
coded  Into  the  state  enumeration  framework  of  the  existing  program. 

An  Inadequate  user  Interface  Is  a  third  limitation  to  PAFT  F77.  An 
interface  enhancement  should  prompt  input  either  from  an  existing  file  or 
manual  CRT  entry.  It  should  have  limited  editing  capabilities  in  order  to 
facilitate  changes  to  existing  Input.  And,  Idiot-proof Ing  should  be  done 
to  check  for  the  correct  data  types  and  sizes  of  descriptive  labels,  prob¬ 
abilities  between  one  and  zero,  and  complete,  connected,  and  acycllcal  fault 
trees . 


i  I  •  y-’fri 
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USER'S  GUIDE 


2.1  Using  PAFT  F77  on  UNIX 

After  having  logged  on  successfully  on  the  UNIX  system,  the  PAFT  F77 
progreua  can  be  run  by  giving  the  command, 

ftree/paft/code  <  {input  file}  >  {output  file) 

ftree/paft/code  specifies  the  location  of  the  PAFT  F77  machine  code  command 
file.  The  first  parameter,  preceded  by  '<'  is  the  file  name  containing 
the  data  that  defines  the  fault  tree  to  be  analysed.  The  second  parameter 
preceded  by  '>'  is  optional.  If  given,  output  from  the  program  is  written 
into  a  file  with  the  given  name.  The  default  is  to  have  output  directed 
back  to  one's  terminal. 

2.2  Input  File  Structure 

The  input  file  structure  to  PAFT  F77  is  a  specific  sequence.  The 
following  list  indicates  this  order,  with  each  number  representing  a 
line  in  the  input  file. 

Fault  tree  name  -  any  alpha-numeric  entry  up  to  60  positions  in  length 
with  no  intermediate  blanks. 

@  Probabilities  or  failure  rates  -  the  character  'p'  or  'f  is 

entered  to  specify  the  input  as  either  in  terms  of  probabilities  or 
failure  rates,  respectively. 

If  failure  rates  are  indicated,  then  lines  @  must  be  entered. 

Q)  Number  of  time  intervals  -  a  digit  is  entered  from  1  to  a  maximum  of  9. 
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@  Time  Intervals  -  real  nonnegative  values  are  entered  separated  by  blanks. 

©  Time  Interval  units  -  alpha-numeric  entry  up  to  8  positions  in  length 
with  no  intermediate  blanks. 

(0  Total  number  of  nodes  and  nutiber  of  replicated  basic  events  -  an  integer 
from  1  to  a  maximum  of  99,  and  an  Integer  from  0  to  a  maximum  of  20, 
respectively,  separated  by  a  blank. 

©  Replicated  node  numbers  (if  any)  -  integer  values  separated  by  blanks. 

Data  @  -  @1  is  entered  for  each  node  in  the  fault  tree,  starting  with  the 

top  node  as  number  1. 

(§)  Node  description,  node  type,  number  of  nodes  above,  number  of  nodes 
below  -  description  is  an  alpha-numeric  entry  up  to  16  positions  in 
length  with  no  intermediate  blanks,  type  is  an  Integer  (see  Table  1 
for  the  index) ,  number  of  nodes  above  and  number  of  nodes  below  are 
Integer  entries  from  0  to  a  maximum  of  20.  All  entries  on  this  line 
are  separated  by  blanks. 

If  the  node  is  a  basic  event,  then  line  must  be  entered. 

Probability  or  failure  rate  -  real  nonnegative  value. 

0  Node  number  adjacencies  -  Integer  entries  of  node  numbers  adjacent  from 
above  then  adjacent  from  below,  separated  by  blanks. 

Repeat  -  0  until  all  nodes  have  been  defined. 

The  following  two  examples  of  input  file  structures  for  the  two  fault 

trees  below  should  prove  helpful. 
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TEST.PAFTCODE 


OR  gate  #1 


AND  gate  #2  AND  gate  #3 


Basic  event  #4  Basic  event  #5  AND  gate  #6 

p  =  0.1  p  *  0.2 


Basic  event  #7  Basic  event  #8 
p«0.5  p*0.6 


i 
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@  Test.paf tcode 
(D  P 
©81 
®  5 

®  top#or  302 
©  2  3 

gatejit2and  2  12 
®  1  A  5 

®  gate//3and  2  12 
<@156 
©  be#A  110 
®  0.1 
®  2 

©  betfS  120 
©  0.2 
<@2  3 

®  gate//6and  2  12 
@378 
©  be//7  110 
©  0.5 
@  6 

©  be//8  110 
©  0.6 
@  6 


INPUT  FILE  FOR:  Test.paf tcode 
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A. compex. fault ■ tree. example 

®  f 
0  3 

®  2.0  3.5  5.0 
hours 
©  22  3 
®  10  11  12 
®  Top  3  0  3 
@  2  3  11 
®  Gate2  213 
@  1  A  10  20 
0  Gate3  312 
@  1  5  11 
0  Gate4  313 
@  2  6  8  12 
0  Gate5  213 
@  3  12  13  14 
0  Gate6  213 
@  4  7  15  16 
0  Gate?  613 
@  6  9  17  18 
0  Gate8  511 
4  10 

0  Gate9  511 
0  7  19 
0  RBEIO  120 
0  0.02 
0  2  8 

0  RBEll  120 
0  0.00004 
0  1  3 

0  RBE12  120 
0  0.0012 


INPUT  FILE  FOR: 


A. compex. fault . tree. example 
(continued) 
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@45 
(D  BE13  110 
@  0.01 
@  5 

®  BE14  110 

(5)  0.9 
@  5 

(D  GatelS  4  12 
0  6  21  22 
d)  BE16  110 
®  0.005 
0  6 

<D  BE17  110 
(9)  0.006 
0  7 

(D  BE18  110 
(9)  0.004 
0  7 

®  BE19  110 

(D  0.01 

©  9 

@  BE20  110 
®  0.001 
©  2 

®  BE21  110 
(D  0.073 
0  15 

®  BE22  110 
®  0.0032 
©  15 


INPUT  FILE  FOR:  A. complex. fault. tree. example 


13 


TABLE  1 

NODE  TYPE  INDEX 


Type  index  // 

1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 


t 


346 

347 


Gate  type 

Basic  event 
AND 
OR 

EXOR 

NOT 

2/3 

x2/3 

2/4 

x2/4 

3/4 

x3/4 

2/5 


a/b 


19/20 

X19/20 


-  5b  +  2a  +  8  {adding  1  for  exclusive  cases 
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2.3  Output  File  Examples 

The  following  two  outputs  are  the  results  of  running  Test.paftcode 
and  A. complex. fault. tree. example. 

Note  that  In  the  second  example  the  marginal  occurrence  rates  of  the 
basic  events  and  the  occurrence  rates  for  the  top  event  are  Invalid  due  to 
the  presence  of  NOT  gate  #8,  NOT  gate  it9  and  EXOR  gate 
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Fault  tree  name:  Test.paftcode 

Total  number  of  events:  8  Number  of  replicated  events:  1 


Node// 

Description 

Type 

Probability 

Above//  /  Below// 

1 

top//or 

OR 

0.74000e-01 

/  2 

3 

2 

gate//2and 

AND 

0.20000e-01 

1/  4 

5 

3 

gate//3and 

AND 

0.60000e-01 

1/  5 

6 

4 

be)/4 

BE 

O.lOOOOe+OO 

2/ 

R- 

5 

be//5 

BE 

0.20000e+00 

2  3/ 

6 

gate//6and 

AND 

0.30000e-H)0 

3/  7 

8 

7 

be// 7 

BE 

0.50000e-H)0 

6/ 

8 

be//8 

BE 

0.60000e+00 

6/ 

Marginal 

Node// 

Description 

Type 

Probability 

Importance 

1 

top/for 

OR 

0.74000e-01 

0 . 10000e401 

R- 

5 

be//5 

BE 

0.20000e+00 

0.37000e+00 

4 

be//4 

BE 

O.lOOOOe+OO 

0.14000e+00 

7 

be// 7 

BE 

0.50000e+00 

0 . 10800e+00 

8 

be//8 

BE 

0.60000e+00 

0.90000e-01 

2 

gate//2and 

AND 

0.20000e-01 

3 

gate//3and 

AND 

0.60000e-01 

6 

gate//6and 

AND 

0.30000e+00 

OUTPUT  FILE  FOR:  Test.paftcode 
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Fault  tree  name:  A. complex. fault. tree. example 

Total  number  of  events:  22  Number  of  replicated  events:  3 

Time  Intervals:  hours  2.00000  3.50000  5.00000 


Node# 

Description 

Type 

Fail  Rate 

Above#  /  Below# 

1 

Top 

OR 

/  2 

3  11 

2 

Gate2 

AND 

1/  4 

10  20 

3 

Gate3 

OR 

1/  5 

11 

4 

Gate4 

OR 

2/  6 

8  12 

5 

Gates 

AND 

3/  12 

13  14 

6 

Gated 

AND 

4/  7 

15  16 

7 

Gate7 

2/3 

6/  9 

17  18 

8 

Gates 

NOT 

4/  10 

9 

Gate9 

NOT 

7/  19 

R-  10 

RBEIO 

BE 

0.20000e-01 

2  8/ 

R-  11 

RBEll 

BE 

0.40000e-04 

1  3/ 

R-  12 

RBE12 

BE 

0.12000e-02 

4  5/ 

13 

BE13 

BE 

O.lOOOOe-01 

5/ 

14 

BE14 

BE 

0.90000e+00 

5/ 

15 

GatelS 

EXOR 

6/  21 

22 

16 

BE16 

BE 

0.50000e-02 

6/ 

17 

BE17 

BE 

0.60000e-02 

7/ 

18 

BE18 

BE 

0.40000e-02 

7/ 

19 

BE19 

BE 

O.lOOOOe-01 

9/ 

20 

BE20 

BE 

O.lOOOOe-02 

2/ 

21 

BE21 

BE 

0.73000e-01 

15/ 

22 

BE22 

BE 

0.32000e-02 

15/ 

Time  Interval: 

2.00000 

hours 

Occurrence  rate 

of  top 

event : 

0.86746e-04 

Marginal 

Marginal 

Node# 

Description 

Type 

Probability 

Importance 

Occur  Rate 

1 

Top 

OR 

0.11980e-03 

O.lOOOOe+01 

0.86746e-04 

R-  11 

RBEll 

BE 

0.79997e-04 

0.99996e-M)0 

0.40000e-04 

R-  12 

RBE12 

BE 

0.23971e-02 

0.16604e-01 

0.19879e-04 

13 

BE13 

BE 

0.19801e-01 

0.20006e-02 

0.19612e-04 

20 

BE20 

BE 

0.19980e-02 

0.93493e-04 

0.93317e-07 

14 

BE14 

BE 

0.83470e+00 

0.47459e-04 

0.70612e-05 

R-  10 

RBEIO 

BE 

0.39211e-01 

0.47640e-05 

0.91555e-07 

17 

BE17 

BE 

0.11928e-01 

0.74371e-06 

0.44095e-08 

18 

BE18 

BE 

0.79681e-02 

0.74080e-06 

0.29399e-08 

16 

BE16 

BE 

0.99502e-02 

0.21311e-06 

0.10551e-08 

21 

BE21 

BE 

0.13584e+00 

0.14901e-07 

0.94014e-09 

22 

BE22 

BE 

0.63796e-02 

0.10994e-07 

0.34961e-10 

19 

BE19 

BE 

0.19801e-0i 

-0.15069e-07 

-0.14773e-09 

2 

Gate2 

AND 

0.18992e-06 

3 

Gate! 

OR 

0.11961e-03 

4 

Gate4 

OR 

0.96088e-K)0 

5 

Gates 

AND 

0.39620e-04 

6 

Gated 

AND 

0.27135e-04 

7 

Gate7 

2/3 

0.19411e-01 

8 

Gates 

NOT 

0.96079e-H)0 

9 

Gate9 

NOT 

0.98020e+00 

15 

GatelS 

EXOR 

0.14049e+00 

OUTPUT  FILE  FOR: 


A. complex .fault . tree . example 
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Time  Interval: 

3.50000 

hours 

Occurrence  rate 

of  top 

event : 

0.12454e-03 

Marginal 

Node// 

Description 

Type 

Probability 

Importance 

1 

Top 

OR 

0.27894e-03 

O.lOOOOe+01 

R-  11 

RBEll 

BE 

0.13999e-03 

0.99986e+00 

R-  12 

RBE12 

BE 

0.4l912e-02 

0.33144e-01 

13 

BE13 

BE 

0.34395e-01 

0.40101e-02 

20 

BE20 

BE 

0.34939e-02 

0.28297e-03 

14 

BE14 

BE 

0.95715e+00 

0.14410e-03 

R-  10 

RBEIO 

BE 

0.67606e-01 

0.14624e-04 

17 

BE17 

BE 

0.20781e-01 

0.63009e-05 

18 

BE18 

BE 

0.13902e-01 

0.62585e-05 

16 

BE16 

BE 

0.17348e-01 

0.18094e-05 

21 

BE21 

BE 

0.22547e+00 

0.13248e-06 

22 

BE22 

BE 

0.11138e-01 

0.74386e-07 

19 

BE19 

BE 

0.34395e-01 

-0.22557e-06 

2 

Gate2 

AND 

0.10214e-05 

3 

Gate  3 

OR 

0.27795e-03 

4 

Gate4 

OR 

0.93269e+00 

5 

Gate5 

AND 

0.13798e-03 

6 

Gate6 

AND 

0.13347e-03 

7 

Gate7 

2/3 

0.33222e-01 

8 

Gate8 

NOT 

0.932396+00 

9 

Gate9 

NOT 

0.96561e+00 

15 

Gatel5 

EXOR 

0.23159e+00 

Time  Interval: 

5.00000 

hours 

Occurrence  rate 

of  top 

event : 

0.15866e-03 

Marginal 

Node# 

Description 

Type 

Probability 

Importance 

1 

Top 

OR 

0.49130e-03 

O.lOOOOe-01 

R-  11 

RBEll 

BE 

0.19998e-03 

0.99971e+00 

R-  12 

RBE12 

BE 

0.59820e-02 

0.48671e-01 

13 

BE13 

BE 

0.48771e-01 

0.59116e-02 

20 

BE20 

BE 

0.49875e-02 

0.57552e-03 

14 

BE14 

BE 

0.98889e+00 

0.29155e-03 

R-  10 

RBEIO 

BE 

0.95163e-01 

0.30163e-04 

17 

BE17 

BE 

0.29554e-01 

0.23738e-04 

18 

BE18 

BE 

0.19801e-01 

0.23514e-04 

16 

BE16 

BE 

0.24690e-01 

0.68309e-05 

21 

BE21 

BE 

0.30580e+00 

0.52343e-06 

22 

BE22 

PE 

0.15873e-01 

0.20998e-06 

19 

BE19 

BE 

0.48771e-01 

-0.12255e-05 

2 

Gate2 

AND 

0.30079e-05 

3 

Gate  3 

OR 

0.48843e-03 

4 

Gate4 

OR 

0.90544e+00 

5 

Gates 

AND 

0.28851e-03 

6 

Gate6 

AND 

0.357563-03 

7 

Gate7 

2/3 

0.46421e-01 

8 

Gate8 

NOT 

0.90484e400 

9 

Gate9 

NOT 

0.95123e+00 

15 

GatelS 

EXOR 

0.31197e+00 

OUTPUT  FILE  FOR:  A. complex. fault . tree . example 


Marginal 

Occur  Rate 

0.12454e-03 

0.40000e-04 

0.39618e-04 

0.38732e-04 

0.28206e-06 

0.55590e-05 

0.27278e-06 

0.37030e-07 

0.24693e-07 

0.88926e-08 

0.74924e-08 

0.23545e-09 

■0.21788e-08 


Marginal 

Occur  Rate 

0.15866e-03 

0.40000e-04 

0.58084e-04 

0.56260e-04 

0.57293e-06 

0.29164e-05 

0.54612e-06 

0.13828e-06 

0.92238e-07 

0.33328e-07 

0.26538e-07 

0.66159e-09 

0.11663e-07 


w 


1 
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PAFT  F77  LISTING  DESCRIPTION 

3.1  Overview  of  Program  Structure 

PAFT  F77  is  comprised  of  fifteen  Fortran  F77  modules.  In  alphabetical 

order  these  are: 

Subroutine  bln(l)  -  adds  1  to  the  binary  array  bln(l)  through  a  recursive 
procedure . 

Function  calpb(n)  -  calculates  and  returns  the  probability  of  node  n  . 

Function  calspb(n)  -  calculates  and  returns  the  probability  of  node  n 
during  state  enumeration. 

Subroutine  Import  -  calculates  the  marginal  Importance  of  all  non-repllcated 
basic  events. 

Subroutine  Input  -  reads  the  Input  file. 

Program  paft  -  the  main  module,  performs  subroutine  calls  and  data  inltillza- 
tion. 

Subroutine  match  -  sets  flags  for  nodes  dependent  on  replicated  basic  events. 

Subroutine  occur (ndone)  -  calculates  the  marginal  occurrence  rates  of  all 
basic  events  and  the  occurrence  rate  of  the  top  event  for  the  ndone^^ 
time  interval. 

Subroutine  output (ndone)  -  outputs  data  for  a  probability  Input,  or  for  the 
ndone^^  time  Interval  of  a  failure  rate  Input. 

Subroutine  qsort(m,n)  -  sorts  the  array  S  from  subscript  m  to  n  ;  S 
Indexes  the  marginal  Importances  of  basic  events  In  decreasing  order. 
Qsort  uses  a  recursive  procedure. 

Subroutine  rduce(l)  -  reduces  the  fault  tree  by  calculating  all  gate  prob¬ 


abilities  not  dependent  on  any  replicated  basic  events.  Rduce  uses 
a  recursive  procedure  starting  with  node  1. 
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Subroutine  remark(l)  -  temporarily  resets  the  flags  signifying  nodal 

dependence  on  replicated  basic  events,  for  the  nodes  above  non-repli- 
cated  basic  event  1  . 

Subroutine  srduce(l)  -  reduces  the  fault  tree  by  calculating  all  gate  prob¬ 
abilities  during  a  specific  state  enumeration.  Srduce  uses  a  recursive 
procedure  starting  with  node  1. 

Subroutine  states  -  calculates  all  gate  prouabilitles  dependent  on  replicated 
basic  events  through  the  state  enumeration  method.  Also  calculates  the 
marginal  Importances  of  all  replicated  basic  events. 

Function  wrdtyp(i)  -  a  character  function  returning  the  alphanumeric  equiva¬ 
lent  of  the  type  index  i  . 

The  diagram  below  shows  the  structure  of  the  procedural  dependence 

between  the  modules. 


wrdtypj  I  occur  |  |  qsort 

\  ♦  / 

[output 


calpb  I 

rduce 


calspb 


( 
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Two  connnon  data  blocks  are  used  to  define  global  variables;  nn  ,  for 
numerical  data,  and  cc  ,  for  character  data. 

3.2  Global  Variable  Definitions 

The  following  list  gives  the  global  variable  definitions  for  PAFT  F77. 
They  are  ordered  according  to  their  position  in  the  common  blocks  cc  and 
nn  . 

Common  block  cc  - 

name  -  the  fault  tree  name;  up  to  60  characters. 
descrp(99)  -  node  descriptions;  up  to  16  characters  each, 
units  -  time  units  for  failure  rate  input;  up  to  8  characters, 
pf  -  probability  or  failure  rate  input  indicator;  1  character. 

Common  block  nn  - 

type(99)  -  node  type  indexes;  integer. 
above(99)  -  number  of  nodes  above;  integer. 
below(99)  -  number  of  nodes  below;  integer. 

nodes(99,20)  -  adjacent  node  numbers  above  and  then  below;  integer. 
prb(99)  -  probability  of  nodes;  real. 

sprb(99)  -  scratch  probability  of  nodes  during  state  enumeration;  real. 

prt(99)  -  marginal  importance  of  nodes;  real. 

fail(99)  -  basic  event  failure  rates;  real. 

time(9)  -  time  intervals;  real. 

sumocr  -  occurrence  rate  of  top  event;  real. 

ocr(99)  -  marginal  occurrence  rates  of  basic  events;  real. 

total  -  total  number  of  nodes;  integer. 
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nrep  -  number  of  replicated  basic  events;  integer, 
nbe  -  number  of  basic  events;  integer. 

replc(20)  -  replicated  basic  event  node  numbers;  integer, 
numti  -  nuiii>er  of  time  Intervals;  integer. 

mark(99)  -  flag  signifying  nodal  dependence  on  a  replicated  basic 
event  or  an  enumerated  basic  event;  integer. 
store(99)  -  scratch  store  for  mark(99) ;  integer 

bn(20)  -  binary  array  signifying  states  of  replicated  basic  events 
Integer. 

dpend(99)  -  flag  signifying  nodal  dependence  on  a  replicated  basic 
event  through  two  or  more  paths;  integer. 
s(99)  -  index  for  order  of  marginal  importances;  integer. 


^  iji 
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APPENDIX 


subroutine  bin(l) 

Integer  bn(20) 

common  /nn/  type(99) ,above(99) ,below(99) ,nodes(99,20) ,prb(99) ,sprb(99) , 
&  prt(99) ,fail(99) ,time(9) ,sumocr,ocr (99) , total, nrep,nbe, 

&  replc(20) ,numti,mark(99) ,store(99) ,bn(20) ,dpend(99) ,s(99) 

c  change  i  th  member  of  bn,  0  to  1  or  1  to  0 

If  (bn(l).eq.O)  then 
bn(i)=l 

c  if  0  to  1  return 

return 

else 

bn(i)=0 

c  if  1  to  0  change  1+1  th  member  of  bn 

call  bin(i+l) 
end  if 
return 
end 

function  calpb(n) 
real  p,pl 

integer  a,b,c,d,fin,start,num,bn(20) ,type(99) ,above(99) ,below(99) , 

&  nodes (99, 20) 

common  /nn/  type(99) ,above(99) ,below(99) ,nodes(99 ,20) ,prb(99) ,sprb(99) , 
prt(99) ,fail(99) ,time(9) ,sumocr,ocr(99) , total, nrep.nbe, 
replc(20) ,numti,mark(99) ,store(99) ,bn(20) ,dpend(99) ,s(99) 

a=above(n) 
b=below(n) 

c  type  of  gate  defines  calpb  operation 

if  (type(n) ,ge. 6)  go  to  50 
p=prb (nodes (n , a+1) ) 
c-2 

go  to  (10,20,30,40)  type(n)-l 
c  AND  gate 

10  do  1  i-a+2,  a+b 

p»p*prb (nodes (n , i) ) 

1  continue 
go  to  70 

c  OR  gate 

20  c*l 

c  EXOR  gate 

30  do  2  j»a+2,  a+b 

pl«prb(nodes(n, j) ) 
p-p+pl-c*p*pl 

2  continue 
go  to  70 

c  NOT  gate 

40  P“l-p 

go  to  70 

c  all  other  gates 

50  p-0.0 
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c  change  type  index  to  (x)a/b  form 

d“(type(n)-b*b+5*b-8) /2 
do  3  g=l,  b 
bn(g)*0 

3  continue 

c  set  up  correct  iteration  bounds  for  generating  function  operations 

if  (mod(type(n) ,2) .eq.l)  then 
start»d 
fin*d 

else 

if  (d.le.b/2)  then 
start=0 
fin=d-l 

else 

started 
fin=b 
end  if 
end  if 

c  generating  function  operation 

do  4  k=l,  b**2 
num=0 

do  5  1=1,  b 

if  (bn(l).eq.l)  num=numfl 

3  continue 

if  ((num.lt.start) .or. (num.gt.fin))  go  to  60 
pl=1.0 

c  probability  calculation 

do  6  in=l,  b 

if  (bn(m).eq.l}  then 
pl*pl*prb  (nodes  (n,a-hn)  ) 
else 

pl=pl* (1-prb (nodes (n,a+m) ) ) 
end  if 
6  continue 

p*p+pl 

60  if  (k.lt.b**2)  call  bin(l) 

4  continue 

if  (d.le.b/2)  p=l-p 
70  calpb=p 

return 
end 

function  calspb(n) 
real  p,pl 

integer  a,b,c,d,fln,start,num,bn(20) ,type(99) ,above(99} ,below(99) , 

&  nodes (99, 20) 

common  /nn/  type(99) ,above(99) ,below(99) ,nodes(99,20) ,prb(99) ,8prb(99) , 
&  prt(99) ,fail(99) ,tlme(9) ,sumocr ,ocr(99) , total, nrep,nbe, 

&  replc(20) ,numti,mark(99) ,store(99) ,bn(20) ,dpend(99) ,s(99) 

a»above(n) 
b=below(n) 
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c  type  of  gate  defines  calspb  operation 

if  (type(n) .ge.6)  go  to  50 
p*sprb (nodes (n , a+1) ) 
c-2 

go  to  (10,20,30,40)  type(n)-l 
c  AND  gate 

10  do  1  i*a+2,  a+b 

p=p*sprb (nodes (n , i) ) 

1  continue 
go  to  70 

c  OR  gate 

20  c=l 

c  EXOR  gate 

30  do  2  j  =a+2 ,  a+b 

pl=sprb (nodes (n , j ) ) 
p=p+pl-c*p*pl 

2  continue 
go  to  70 

c  NOT  gate 

40  P=1“P 

go  to  70 

c  all  other  gates 

50  p-0.0 

c  change  type  index  to  (x)a/b  form 

d=(type(n) -b*b+5*b-8) /2 
do  3  g=l,  b 
bn(g)=0 

3  continue 

c  set  up  correct  iteration  bounds  for  generating  function  operations 

if  (mod(type(n) ,2) .eq.l)  then 
start=d 
fin=d 

else 

if  (d.le,b/2)  then 
start=0 
fin-d-1 

else 

start-d 
fin»b 
end  if 

c  generating  function  operation 

end  if 

do  4  k>=l,  b**2 
num-0 

do  5  1-1,  b 

if  (bn(l).eq.l)  num-num+l 
5  continue 

if  ((num.lt.start) .or. (num.gt.fln))  go  to  60 
pl-1.0 

c  sprob  calculations 

do  6  m-1,  b 
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if  (bn(m).eq.l)  then 

pl=pl*sprb (nodes (n,a+m) ) 
else 

pl=pl* ( 1-sprb (nodes (n , a+m) ) ) 
end  if 
continue 

p=p+pl 

If  (k.lt.b**2)  call  bln(l) 
continue 

if  (d.le.b/2)  p“l-p 

calspb=p 

return 

end 

subroutine  Import 

integer  l,c,type(99) , above (99) .total, nrep,replc( 20) ,bn(20) , 
store(99) .mark (99) 

real  p.x 

common  /nn/  type (99) .above (99) .below (99) .nodes (99 ,20) ,prb(99) ,sprb(99) , 
prt(99) ,fall(99) ,time(9) .sumocr ,ocr(99) , total, nrep.nbe, 
replc(20) ,numti,mark(99) ,store(99) ,bn(20) ,dpend(99) ,s(99) 
top  event  always  has  a  marginal  Importance  of  1.0 
prt(l)=1.0 

reclaim  original  mark  flags 
do  1  i=l,  total 
store(i)*mark(i) 
continue 

select  only  nonreplicated  basic  events 
do  2  n=l,  total 

if  (type(n) .gt.l)  go  to  2 
if  ( above (n) .gt.l)  go  to  2 
call  remark (n) 
sprb(n) =0.0 

basic  event  first  off  then  on 
do  3  m=l,  2 

if  (m.eq.2)  sprb(n)=1.0 
do  A  g=l,  nrep 
bn(g) =0 
continue 
c=2**nrep 

state  enumeration  on  nonreplicated  basic  event 
do  5  j“l,c 

p-1.0 

do  6  k=l,  nrep 
l«replc(k) 
sprb(l)=bn(k) 
probability  calculation 

if  (bn(k).eq.l)  then 
p=p*prb(l) 

else 

p.p*(l_prb(l)) 
end  if 
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6  continue 

call  srduce(l) 
x«p*sprb(l) 

c  marginal  Importance  calculation 

If  (m.eq.2)  then 
prt (n) “prt (n)+x 
else 

prt(n)*prt(n)-x 
end  If 

If  (j.lt.c)  call  bln(l) 

5  continue 

3  continue 

do  7  h=l,  total 

If  (store(h) .eq.l)  sprb(h)=prb(h) 

7  continue 

2  continue 

return 

end 

subroutine  Input 

Integer  type(99) ,above(99) ,below(99) ,nodes(99 ,20) , total, numtl, 

6  nrep,replc(20) 

character  name*60 , descrp (99) *16 ,unlts*8 ,pf 

common  /nn/  type(99) ,above(99) ,below(99) ,nodes(99 ,20) ,prb(99) ,sprb(99) , 
&  prt(99) ,fall(99) ,tlme(9) ,sumocr,ocr (99) , total, nrep,nbe, 

&  replc(20) ,numtl,mark(99) ,store(99) ,bn(20) ,dpend(99) ,s(99) 

common  /cc/  name, descrp(99) , units, pf 
read (5,*)  name 
read (5,*)  pf 
If  (pf.eq.’p*)  go  to  10 
read (5,*)  numtl 

read(5,*)  (tlme(k) ,k«l,  numtl) 
read(5,*)  units 
10  read (5,*)  total, nrep 

If  (nrep.gt.O)  read(5,*)  (replc(n) ,n“l,  nrep) 
do  1  1~1,  total 

read(5,*)  descrp(l) ,type(i) ,above(l) ,below(l) 

If  (type(l) .eq.l)  then 
read(5,*)  prb(l) 

c  fall  to  be  changed  to  a  probability  In  main  routine 

fall(l)-prb(l) 
end  If 

read(5,*)  (nodes (1,J) ,j“l,  above(l)+below(l) ) 

1  continue 

return 
end 

program  paft 

Integer  j ,n,type(99) .total, nrep, above(99) ,below(99) ,replc(20) , numtl, 

&  nbe,nodes(99,20),mark(99),B(99) ,bn(20) ,dpend(99) ,store(99) 

real  prb(99) ,sprb(99) ,prt(99) ,fall(99) ,tlme(9) .sumocr , 

&  ocr(99) 
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character  nanie*60 ,  descrp  (99)  *16 , units*8  ,pf 

common  /nn/  type (99) .above (99) ,below(99) .nodes (99, 20) ,prb(99) ,sprb(99) , 
&  prt(99) ,fail(99) ,time(9) ,sumocr,ocr(99) , total, nrep.nbe, 

&  replc(20)  ,numtl.mark.(99)  ,store(99)  ,bn(20)  ,dpend(99)  ,s(99') 

common  /cc/  name, descrp(99) .units, pf 
call  Input 

c  prt(top)  Is  always  equal  to  1.0 

s(l)*l 
nbe=l 

c  j  keeps  count  of  non-basic  events 

j=total 

c  set  up  s  array  for  correct  Input  to  qsort  and  output  routines 

do  1  k=l,  total 
dpend(k)*0 
mark(k)«l 

if  (k.eq.l)  go  to  1 
n*total+2-k 
if  (type(n) .eq.l)  then 
nbe=nbe+l 
s(nbe)*n 

else 

s(j)=n 

j=j-l 

end  if 

1  continue 
call  match 

c  number  of  output  loops 

if  (pf.eq.'p')  numti=l 
do  2  ndone-1,  numtl 

c  Initlllzatlons;  fall  into  prb.prb  into  sprb.and  prt 

do  3  i»l,  total 

if  ( type (i) .eq.l)  then 

if  (pf.eq.'f’)  prb(l)*l-exp(-fail(i)*time(ndone)) 
else 

prb(i)=0.0 
end  if 

8prb(i)»prb(i) 

prt(l)-0 

3  continue 

call  rduce(l) 
if  (nrep.eq.O)  go  to  10 
call  states 
10  call  import 

call  output (ndone) 

2  continue 
stop 
end 

subroutine  match 

integer  l,c,n, J ,llsc(99) ,above(99) ,nodes(99,20) ,dpend(99) ,nrep, 
mark (99) ,replc(20) 


& 


CSUiH  (SU<H  b9«<0 
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common  /nn/  type(99) ,above(99) ,below(99) ,node8(99,20) ,prb(99) ,sprb(99) , 
prt(99) ,fall(99) ,time(9) ,sumocr,ocr(99) , total, nrep,nbe, 
replc(20) ,numti,mark(99) ,store(99) ,bn(20) ,dpend(99) ,s(99) 
do  1  i*l,  nrep 

list  stores  replicated  basic  events  already  processed 
list(l)=replc(i) 

1=1 

c=l 

0  do  2  k=l,  above (list(c)) 

process  starts  at  replicated  basic  event  and  marks  up  the  fault  tree 
n=nodes(list(c) ,k) 

j=l 

0  if (llst(j) .eq.n)  dpend(n)=l 

if  (dpend(n) .eq.l)  go  to  2 

j=j+l 

if  (l.le.l)  go  to  10 
1=1+1 
list(l)=n 
mark(n)=0 
continue 

next  replicated  basic  event 
c=c+l 

if  (c.le.l)  go  to  20 
continue 
return 
end 

subroutine  occur (ndone) 

Integer  type(99) , total 
real  q 

common  /nn/  type(99) ,above(99) ,below(99) ,nodes(99,20) ,prb(99) ,sprb(99)  , 
prt(99) ,fall(99) ,time(9) ,sumocr,ocr(99) .total, nrep, nbe, 
replc(20) ,numti,mark(99) ,8tore(99) ,bn(20) ,dpend(99) ,s(99) 

sumocr*0.0 
q-l-prb(l) 
do  1  k-2,  total 

if  ( type (k) .eq.l)  then 

ocr (k) « ( 1-prb (k) ) *  f  ail (k) *pr t (k) /q 
sumocr=sumocr+ocr (k) 
end  if 
1  continue 

ocr (l)-8umocr 

return 

end 

subroutine  output (ndone) 

Integer  a,b,type(99) ,above(99) ,below(99) , node8(99, 20) .total, numti, 

&  nrep, s(99) , nbe 

character  name*60 , descrp (99) *16 ,b lanks*!! , uni ts*8,word*6 ,wrdtyp*6 ,r*2 ,pf 
common  /nn/  type(99) ,above(99) ,below(99) ,node8(99 ,20) ,prb(99) ,8prb(99) , 

&  prt(99) ,fall(99) ,time(9) ,sumocr,ocr(99) .total, nrep, nbe, 

6  replc(20) , numti, mark(99) ,store(99) ,bn(20) ,dpend(99) ,s(99) 
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common  /cc/  name, descrp (99) .units, pf 
blanks*' 

c  go  to  time  Interval  output  code  if  ndone  Is  greater  than  1 

if  (ndone. gt.l)  go  to  1000 
open(6, form*' print') 
write (6, 10)  name 

10  format ( 'OFault  tree  name:  ',a60) 

write (6, 20)  total. nrep 

20  formate  Total  number  of  events:  ',12,'  Number  of  replicated', 

&  '  events:  ',12) 

c  go  to  failure  rate  output  code  if  input  is  in  terms  of  failure  rates 

if  (pf.eq.'f)  go  to  2000 
c  probability's  initial  output 

write(6,30) 

30  format ( 'ONode#  Description  Type  Probability  ', 

&  'Above//  /  Below//') 

do  1  k=l,  total 
a*above (k) 
b=below(k) 
r=' 

if  (a. gt.l)  r*'R-' 

word*wrdtyp ( type (k) ) 

write (6 , 40)  r , k, descrp (k) .word , prb (k) 

40  format(lx,a2,i3,3x,al6,lx,a4,2x,el2.5,$) 

if  (a.eq.O)  wrlte(6,50)  blanks 
if  (a.eq.l)  wrlte(6,60)  nodes(k,l) 
if  (a.eq.2)  write(6,70)  nodes(k,l),nodes(k,2) 
if  (a.qt.2)  write(6,80)  (nodes (k,i) ,1*1,  a) 

30  format(lx,all,$) 

60  format (9x, 13, $) 

70  format (6x, 213,$) 

80  format(3x, 20(13,$)) 

write (6, 90)  (nodes (k,j) ,j»a+l,  a+b) 

90  format(lx, ' / ' ,20i3) 

1  continue 

c  probability's  secondary  output 

write (6, 100) 

100  format ('0  Marginal') 

write(6,110) 

110  format ('  Node//  Description  Type  Probability  Importance') 

call  qsort(2,nbe) 
do  2  n*l,  total 
r*' 

if  ( above (s(n)). gt.l)  r*'R-' 
word*wrd typ (type(s(n))) 
if  ((n,eq.l) .or. (type(8(n)) .eq.l))  then 
wrlte(6,120)  r,s(n) ,descrp(8(n)) ,word,prb(8(n)) ,prt(8(n)) 

el8e 

write (6 , 120)  r , 8 (n) , de8crp (s (n) ) .word , prb (8 (n) ) 
end  if 

format(lx,a2,13,3x,al6,lx,a4,2x,el2.5,2x,el2.5,2x,el2.5) 
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2  continue 
return 

c  failure  rate's  Initial  output 

2000  wrlte(6,130)  units, (tlme(f) ,f«l,  numti) 

130  formate  Time  Intervals:  ',a8,3f8.5) 

write (6,140) 

140  f ormat ( ' ONode//  Description  Type  Fail  Rate  Above//', 

&  '  /  below#') 

do  3  1=1,  total 
a*above(l) 
b=below(l) 
r*' 

if  (a.gt.l)  r='R-' 
word-wrdtyp(type(l) ) 
if  (type(l) .eq.l)  then 
write (6 , 150)  r , 1, descrp (1) .word , fail (1) 
else 

wrlte(6,160)  r,l,descrp(l) .word, ' 
end  if 

150  format(lx,a2,13,3x,al6,lx,a4,2x,el2,5,$) 

160  f  ormat ( lx , a2 , 13 , 3x, al6 , lx , a4 , al4 , $ ) 

if  (a.eq.O)  wrlte(6,50)  blanks 

if  (a. eq.l)  wrlte(6,60)  nodes(l,l) 

if  (a.eq.2)  write(6,70)  nodes(l,l) ,nodes(l,2) 

if  (a.gt.2)  write(6,80)  (nodes (l.m) ,m=l,  a) 

write(6,90)  (nodes (l,g) ,g*«a+l,  a+b) 

3  continue 

c  failure  rate's  secondary  output 

1000  write (6, 170)  time (ndone) .units 

170  formatCOTlme  interval:  '  ,f8.5,lx,a8) 

call  occur (ndone) 
write (6, 180)  sumocr 

180  formatC  occurrence  rate  of  top  event:  ',el2.5) 
write(6,190) 

190  formate  Marginal  ', 

&  '  Marginal ' ) 

write (6, 200) 

200  formatC  Node#  Description  Type  Probability  Importance', 

&  '  Occur  Rate') 

call  q8ort(2,nbe) 
do  4  h=l,  total 
r=' 

if  (above (s (h) ) .gt.l)  r='R-' 
word=wrdtyp(type(8(h))) 
if  ((h. eq.l). or. (type(8(h)). eq.l))  then 
wrlte(6,120)  r,8(h) ,de8crp(8(h)) , word, prb (8(h)) ,prt(8(h)) ,ocr(8(h)) 

else 

write(6,120)  r,8(h) ,descrp(s(h) ) ,word,prb(s(h) ) 
end  if 

4  continue 
return 
end 
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subroutine  qsort(m,n) 

integer  1(99) ,11(99) ,i(99) , j j (99) ,L,k,s(99) .total 
real  x 

connnon  /nn/  type(99) ,above(99) ,below(99) .nodes(99,20) ,rrb(99) ,sprb(99) 

prt(99) ,fall(99) ,time(9) ,suinocr,ocr(99) .total, nrep.nbe, 

//  ,/®P^‘'<20),nuinti,mark(99),store(99),bn(20),dpend(99),s(99) 

If  ((in.eq.2).and.(n.eq.total))  1-0 
1  is  recursive  level  indicator 
1-1+1 
ii(l)-m 
jj(l)-n 
KD-m 
j(l)»n 

x=prt(s(j(l))) 

standard  recursive  quicksort  routine 
if  (prt(s(i(l))) .le.x)  go  to  20 
i(l)-i(l)+l 
go  to  10 

if  (prt(s(j(l))).ge.x)  go  to  30 
j(l)-j(l)-l 
go  to  20 
30  if  (l(l).le.j(l))  then 

k-s(l(l)) 
s(l(l))-s(j(l)) 
s(j(l))-k 
l(l)-i(l)+l 
j(l)-J(l)-l 
end  if 

if  (i(l).le.j(l))  go  to  10 
c  resort  all  smaller  than  x 

if  (ii(l).lt.J(l))  then 
call  qsort(ii(l),j(l)) 

1-1-1 
end  if 

c  resort  all  larger  than  x 

if  (i(l).lt.jj(l))  then 
call  qsort(i(l),jj(l)) 

1-1-1 
end  if 
return 
end 


subroutine  rduce(i) 

Integer  "(99) ,k(99) ,l,above(99) .below(99) .nodes(99,20) ,dpend(99) 
common  /nn/  type(99) .above(99) ,below(99) ,nodes(99 ,20) ,prb(99) ,sprb(99), 
prt(99)  ,fail(99)  ,tiiiie(9)  ,8umocr,ocr(99)  , total, nrep.nbe, 

.r  ^«Pl‘^<20),numtl,mark(99),8tore(99),bn(20),dpend(99),s(99) 

if  (l.eq.l)  then 
1-0 

n(l)-l 
end  if 

1  is  recursive  level  Indicator 

1-1+1 
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c  k  Is  Iterative  breadth  pointer 

k(l)«above(n(l))+l 
10  n(l+l)»nodes(n(l) ,k(l)) 

c  recurse  depth  first 

If  (below(n(l+l)) .gt.O)  call  rduce(n(l+l)) 

Iterate  breadth  second 
If  (k(l) .lt.above(n(l))+below(n(l)))  then 
k(l)-k(l)+l 
to  to  10 
end  If 

c  probability  calculations 

If  (dpend(n(l)) .eq.O)  then 
prb(n(l))-calpb(n(l)) 
sprb (n (1) ) »prb (n(l) ) 
end  If 
1-1-1 
return 
end 

subroutine  remark (1) 

Integer  j ,n,store(99) ,nodes(99,20) , mark (99) , total 

common  /nn/  type(99) ,above(99) ,below(99) , nodes (99, 20) ,prb(99) ,sprb(99) , 
&  prt(99) ,fail(99) ,tlme(9) .sumocr ,ocr(99) , total, nrep.nbe, 

&  replc(20) ,numti,mark(99) ,store(99) ,bn(99) ,dpend(99) ,s(99) 

c  store  original  mark  in  store  array 

do  1  k*l,  total 
mark (k) -store (k) 

1  continue 

j-i 

10  n-nodes ( j , 1) 

c  unmark  all  nodes  above  1 

mark(n)-0 
j-n 

if  (j.gt.l)  go  to  10 
return 
end 

subroutine  srduce(l) 

Integer  n(99),k(99),l,above(99) ,below(99) ,nodes(99,20) ,mark(99) 
common  /nn/  type(99) ,above(99) ,below(99) ,node8(99 ,20) ,prb (99) ,8prb(99) , 
prt (99) , fail (99) , time (9) , sumocr , ocr (99 ) , total , nrep ,nbe , 
replc(20) ,numti,mark(99) ,8tore(99) ,bn(20) ,dpend(99) ,8(99) 

If  (1 .  eq .  1)  then 
1-0 

n(l)-l 
end  If 

c  1  is  recursive  level  indicator 
1-1+1 

c  k  Is  Iterative  breadth  pointer 
k(l)-above(n(l))+l 
10  n(l+l)-nodes(n(l) ,k(l)) 

c  recurse  depth  first 

If  (mark(n(l+l)) .eq.O)  call  srduce(n(l+l)) 
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c  Iterate  breadth  second 

If  (k(l) .lt.above(n(l))+below(n(l)))  then 
k(l)=k(l)+l 
go  to  10 
end  if 

c  sprob  calculations 

sprb(n(l))=calspb(n(l)) 

1=1-1 
return 
end 

subroutine  states 

Integer  l,c, type(99) , total, nrep,replc(20) ,bn(20) ,dpend(99) 
real  p,x 

common  /nn/  type(99) ,above(99) ,below(99) ,nodes(99,20) ,prb(99) ,sprb(99)  , 
prt(99) ,fall(99) ,tlme(9) ,sumocr,ocr(99) , total, nrep,nbe, 
replc(20) ,numtl,mark(99) ,store(99) ,bn(20) ,dpend(99) ,s(99) 
do  1  1=1,  nrep 
bn(l)=0 
continue 

state  enumeration  process 
c=2**nrep 
do  2  j=l,  c 

p=1.0 

do  3  k=l,  nrep 
l=replc(k) 
sprb(l)=bnCk) 
c  probability  calculations 

If  (bn(k).eq.l)  then 
p=p*prb(l) 

else 

p=p*(l-prb(l)) 
end  If 

3  continue 

call  srduce(l) 
do  4  n=l,  total 

If  (dpend(n) .eq.O)  go  to  4 
x=p*sprb(n) 
prb(  ^-prb(n)+x 
if  (n.gt.l)  go  to  4 

c  replicated  basic  event  marginal  Importances  calculations 

do  5  m=l,  nrep 
l=replc (m) 

If  (bn(m).eq.l)  then 
prt(l)-prt(l)+x/prb(l) 
else 

prt (1) -prt (l)-x/ (1-^Tb (1) ) 
end  if 

5  continue 

4  continue 

If  (J.lt.c)  call  bln(l) 

2  continue 

return 
end 
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character*6  function  wrdtyp(l) 
integer  a.b.ind 

character  word(5)*6,  cha*2,  chb*2 
data  word/ 'BE' , 'AND' , 'OR' , 'EXOR' , 'NOT'/ 
if  type  index  is  less  than  6  return  correct  data  word 
if  (l.lt.6)  then 
wrdtyp=word (i) 
return 
end  if 

change  type  index  to  (a)a/b  alpha-numeric  form 
b=(sqrt(4*i-23.)+5)/2 
a=(l-b*b+5*b-8)/2 
write(cha, ' (12) ')  a 
write(chb, ' (12) ')  b 
if  (mod(i,2) .ne.O)  then 
•  rdtyp='X' 
ind=2 

else 

ind=l 
end  if 

if  (a.ge.lO)  then 
wrdtyp(lnd: )=cha// ' / ' 
ind=ind+3 

else 

wrdtyp(lnd:)=cha(2:2)//'/' 
ind=ind+2 
end  if 

if  (b.ge.lO)  then 
wrdtyp(ind:)=chb 
ind=ind+2 

else 

wrdtyp(ind:)=chb(2:2) 
lnd=ind+l 
end  if 
return 
end 
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